Synchronization in the BCS Pairing Dynamics as a Critical Phenomenon 
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Fermi gas with time-dependent pairing interaction hosts several different dynamical states. Cou- 
pling between the collective BCS pairing mode and individual Cooper pair states can make the latter 
either synchronize or dephase. We describe transition from phase-locked undamped oscillations to 
Landau-damped dephased oscillations in the collisionless, dissipationless regime as a function of 
coupling strength. In the dephased regime, we find a second transition at which the long-time 
asymptotic pairing amplitude vanishes. Using a combination of numerical and analytical methods 
we establish a continuous (type II) character of both transitions. 



Recent discovery of BCS pairing in fermionic vapors Q, 
2], made possible by control of interactions in trapped 
cold gases |3j, has renewed interest in quantum collective 
phenomena jj]. Advanced detection techniques and long 
coherence times in vapors enable time-resolved studies of 
new collective modes, such as spin waves and the BCS 
pairing mode 

Interaction between a collective mode and constituting 
particles is key for our understanding of dynamics in var- 
ious systems, from plasma to quantum gases. One of the 
most surprising of these phenomena is Landau damping 
which occurs in a collisionless regime via direct dissipa- 
tionless energy transfer from the collective mode to single 
particles. Its nondissipative and thus reversible charac- 
ter leads to a variety of regimes, notably to quench- 
ing of the damping, first explored in plasma physics ||. 
Remarkably, a linearly damped mode can regrow and 
transform to a stationary oscillatory Bernstein-Greene- 
Kruskal mode. This fascinating prediction was confirmed 
experimentally only recently 

Naturally, the richness of these nonlinear phenomena 
makes it tempting to look for their analogs in cold gases. 
Collisionless damping in cold gases was considered, in the 
linear regime, for optical excitations lid . spin waves [ill 
ll2T ] , and excitations in optical lattices yll. Motivated by 
the work on fermion superfluidity El M 15, here we fo- 
cus on the pairing dynamics of fermions |l4l EH. EH. Et| 
induced by a sudden change of interaction. The colli- 
sionless regime becomes practical in this case due to long 
elastic collision times r e i ^> ta = fi/A E3, where A is 
the BCS gap. The pairing mode of a small amplitude 
oscillates at a frequency 2A/H and exhibits collisionless 
dephasing 18]. These conclusions were extended recently 
to the nonlinear regime [l9| ]. 

This behavior changes drastically as the perturbation 
increases. The main result of this work, as summarized 
in Fig^ is prediction of a dynamical transition resulting 
from competition between synchronization and collision- 
less dephasing, taking place as a function of the initial 
pairing gap, A s . We found three qualitatively differ- 
ent regimes (A, B, and C) with the critical points at 
A AB = e _7r/2 A and A B c = e^Ao, where A is the 
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FIG. 1: Three regimes of the pairing dynamics vs. the initial 
gap value A s . In synchronized phase (A), A 3 < Aab, the 
pairing amplitude oscillates between A_ and A+. In the de- 
phased regime (B, C), the pairing amplitude saturates to a 
constant value, A n , when Aab < A s < Abc, and decreases 
to zero at A s > Abc- Dashed line: The stationary gap value 
A(T») reached in a closed system after equilibration. 



equilibrium pairing amplitude in the final BCS state. Be- 
low the A-B transition, A s < Aab, individual Cooper 
pair states synchronize and the pairing amplitude oscil- 
lates between A + and A_ without damping. In con- 
trast, in the interval Aab < A s < Abc the pairing 
amplitude is Landau-damped and exhibits decaying os- 
cilation, saturating at an asymptotic value, A , with non- 
monotonic dependence on A s . A second transition occurs 
at A s > Abc, where the dynamics becomes overdamped, 
and A(t) decreases to zero without oscillations. The os- 
cillation amplitude and the asymptotic value A a vanish 
continuously at the critical points A-B and B-C, as in 
a type II transition. We demonstrate that these results 
are consistent with the spectral analysis [2(| based on the 
integrability of the problem. 

We also address the behavior on a long time scale, 
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t ^ T e i, after dissipation sets in. We find that energy re- 
laxation in a closed system, such as an atom trap, makes 
it evolve to a new equilibrium state. Both the tempera- 
ture T* and the gap A(T») exhibit a nonmonotonic de- 
pendence on the initial conditions (Fig.l). 

In our analysis of the BCS problem we employ the 
well known pseudospin formulation [2l| in which spin 1/2 



operators si? 



describe Cooper pairs (p, — p). 



The BCS Hamiltonian takes the form 



w = -5> p «;-A(*)x; 



(1) 



where e p = p 2 /2m — \x is the free particle spectrum with 
p, the Fermi energy. Here we consider the time evolution 
induced by an instantaneous change of interaction from 
A s at t < to A at t > 0. In the spin formulation, Eq.JTJ, 
the dynamics is of a Bloch form 



dt 



= 2b p x r p , 



bp = -(A x , A„,e p ), 



(2) 



where r p = 2(s p ) are classical vectors, and the effective 
magnetic field b p depends on the pairing amplitude A. 
The latter is defined self-consistently: 



A = A, 



- iA„ 



9 Z^'P- 



(3) 



We first present numerical results for the dynamics J2J, 
(J3J . The Runge-Kutta method of the 4th order was used 
with N = 10 4 , 10 5 equally spaced discrete energy states 
within a band W — 50Ao with a constant density of 
states v{Ep). As an initial state we take 



r+(0) = 



rj(0) 



(4) 



which describes the T = paired ground state |21| . 
Without loss of generality we set A(t) = A x , since the 
phase of A is a constant of motion due to the particle-hole 
symmetry of the model. The interaction constants A s and 
A define, via the self-consistency relation ©, the initial 
and final equilibrium BCS gap values, A s = We~ 1 ^ g % 
g s = \ s v < 1, Ao = We~ x l Q , g = \v, which we use to 
parameterize the system. 

We observe three qualitatively different dynamical 
regimes. The initial states with a relatively small gap 
give rise to undamped oscillations (FigjSJi). In this case 
A(t) oscillates non-harmonically between A_ and A + 
(the regime A in Fig^J. Synchronization of different 
Cooper pair states results from their interaction with the 
mode singled out by BCS instability of the initial state, 
similar to the evolution from the normal state [lil ]. 

Desynchronization takes place at A s > Aab = 0.2lAo 
giving rise to two different regimes exhibiting dephasing, 
underdamped and overdamped (B and C, Fig. ^| . The 
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FIG. 2: (a): The pairing amplitude A(t) for the initial 
state IgJ with A s = 0.05Ao as recorded from the simulation, 
oscillating between A+ = 0.97Ao and A_ = 0.31Aq; Syn- 
chronization (b,c): the phase <f> p time dependence, Eq.Q, for 
e p = 0, Ao,2Aq, and frequency ui p vs. e p , Eq.©. 

former, illustrated in FigOJi, is simplest to understand 
for a small initial deviation, A s ~ Ao by lineariz- 
ing Bloch equations about the equilibrium state. The 
analysis predicts damped oscillations at long times: 

A(i) = A a + A(t) sin(2A a t + a), A{t)(xt~ 1/2 . (5) 

The power-law decay of A(t) was explained in Ref. 0] 
by interaction of the collective mode with the continu- 
ous spectrum of excitations with energies above 2A Q and 
linked to the linear Landau damping. In the spin formu- 
lation, the dephasing results from the Larmor frequency 
of spin precession b p being a continuous function of e p . 
An extension of this argument to the nonlinear regime 
was proposed recently in Ref. 19] which, however, did not 
clarify the range of its validity. The dephased time evolu- 
tion similar to JHJ was also reported in Refs. [lfil ItH . l22j | . 

In the overdamped regime, A s > Abc = 4.8lAo 
(FigJSp), A(t) decays to zero without oscillations. This 
behavior can be understood in the limit A s /Ao ^> 1, 
i.e. when the coupling in the initial paired state in Q is 
suddenly completely turned off. For different spins pro- 
cessing freely and independently one obtains 



r+(0), A(t » A; 1 ) oc (A.t)- 1 /a e -aA.t < 



The fast dephasing can also be understood by noting that 
the energy distribution in (gj corresponds to an effective 
temperature T ~ A s which exceeds T c for Ao (see below) . 

To fully exhibit phase locking in the synchronized 
regime which abruptly disappears in the dephased 
regime, we now explore the phase dynamics. It is con- 
venient to measure precession angles relative to time- 
independent b p = — (A a ,0, e p ), where A Q is the asymp- 
tote A(t — > oo) in the regimes B, C, and the average 
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FIG. 3: Dephased dynamics (a): Simulated A(t) for the initial 
states Q with A s = 0.21Ao, 4.5Ao with the asymptotic val- 
ues A a = O.8IA0, 0.12Ao; Overdamped dynamics (b): Same 
as in (a) with A s = 4.81Aq and A a = 0; (c): The phase <f> p 
time dependence, Eq.@, for energies e p — 0,Ao,2Ao (bot- 
tom to top) for A s = 0.21Ao; (d): The frequency u) p vs. e p 
for A a = 0.21Ao (dashed line) and A s = 4.81A (solid line). 

value of oscillating A(i) in A (dash-dotted line in Fig.^l. 
The angle and frequency of precession are defined by 



-t- X 

rip = n P 



!^ W , u p (t) = d4> p /dt, (6) 



where the vectors n p are obtained from r p by a rota- 
tion about the y axis which maps z onto b p : = r®, 



n p + %n p 



e wp( r x _|_ j rp ) j with the rotation angle de- 
fined by tan# p = A a /|e p |. The phase evolution, which 
becomes linear at long times r 3> A^" 1 (Figs l2l3|) . can be 
characterized by average frequency (phase slope) uj p = 
(d<p p /dt) = (0 p (t)— 0p(O))/r. While in the regime A dif- 
ferent p states phase lock (Fig|2p,c), in the regimes B, C 
the frequencies u> p have dispersion (Fig|3t,d). The latter 
reproduces quasiparticle spectrum, uj p = 2(e p + A 2 ) 1 / 2 
with the long-time asymptote A a which vanishes in the 
overdamped regime C. 

We observe a qualitative change in behavior, with u> p 
dispersing in the regions B 7 C and phase locking in the 
region A. However, the oscillation amplitude ^(A + — A_) 
in A and the asymptotic amplitude A a in B vanish con- 
tinuously at the critical points A = Aab,^bc (FigQJ, 
indicating a type II transition. 

Until now, we considered the dissipationless dynamics 
at times shorter than the quasiparticle relaxation time, 
t ;$ t c i. Using the energy balance argument, one can 
determine the system state at long times, t ^> T e \. To ac- 
count for system equilibration, one needs to consider the 
full many-body Hamiltonian which enables elastic scat- 
tering of individual quasiparticles, omitted in Eq.lfT)l. For 
a closed system, such as an atomic trap, the final tem- 
perature Ts,, and the gap A* can be determined from the 
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FIG. 4: The pairing amplitude A(t) as a function of time 
for the initial state (J1J with A s = O.OlAo as recorded from 
numerical simulation. At long time, numerical A(i) matches 
the analytic form in Eg. 1131 . oscillating non- harmonically be- 
tween A+ = 0.997A and A_ = 0.093A . 

energy conservation condition E t -%, Ttsl (T*) — Eq, where 
Eq is the energy immediately after interaction switching, 



En 



A 2 



A 

2 -T, 



(7) 



with the spectrum e p = (e 2 - 
the energy of the final state: 



A 2 ) 1 / 2 , an d £ t>Tel (T*) is 



E, 



Here n p 
T = T* 



e P (A*) 



E 
p 



[e p - (1 -2n p )e p (A»)] 



A| 
A ' 



(8) 



e e p (A,)/IV 

= (4 + a 2 : 



describes equilibrium with 
1 ' 2 . After integrating over 



e p , we arrive at the equation for X 1 ,: 



F 



A* 
2Tl 



= 1 - 



In 



A. 
A 



(9) 



where F(u) = 2 L°° dxcosh.2x [1 — tanh (ucoshx)], and 
a = 1 — gln(A s /A ). From Eq.© we obtain T* and 
the equilibrium gap A* = A(T»). They are approxi- 
mately constant, with T» « 0.72T C , A(T») « 0.81A in 
the regime A, are described by a non-monotonic function 
in B, and vanish in C . The system turns normal in the 
final state for A s > f(g)A , f(g) = 2.2 + 0.86.g + 0{g 2 ). 
Numerical solution of Eq.© at g = 0.26 is displayed in 
Fig^ (dashed line). 

It is instructive to compare our results to the spectral 
analysis based on the integrability of the BCS Hamilto- 
nian 20, 23, 24] . There is an infinite number of commut- 
ing integrals of motion, R n = L p s p , parameterized by p, 
where following Ref. |2fJ we employ the Lax vector, 



L p = z + A ^2 ■ 
p'#p 



(10) 
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We will need L 2 = 4(L p s p ) 2 which is also conserved (due 
to the Pauli matrices algebra combined with Ea. (|10[l h 
The mean-field expressions are obtained by substituting 
the averages (s p ) instead of s p . 

The spectral polynomial defining the evolution of in- 
dividual states is proportional to the square of the Lax 
vector |1>], Q(e) = L 2 (e) U p (^-^p) 2 , where L(e p ) = L p . 
The pairs of complex roots of the spectral equation 
Q(e) — uniquely determine the long-time dynamics 
of the system 20]. Evaluation of L 2 (y) for the initial 
state (@J is straightforward: 



We are grateful to Boris Spivak and Kumar Raman for 
helpful discussions. 



ln^ + yG(y) 



(11) 



where y = e/A s and 



1 f c 



dx 



y-sinhx 2y/l + y 2 "y- ^l + y 2 ' 



In- 



where y is complex. One can show that the roots of L 2 (y) 
lie on the imaginary axis. Upon variation of A s / Aq they 
disappear at y = 0. Expanding about y = 0, we obtain 



V{y)/g 2 



ki- 



ln 



As 
A 



0(y), (12) 



Thus L 2 (y) has no complex roots at A s /Ao > e 7 ^ 2 , one 
pair of roots y = ±iu when e^ 7r/ ' 2 < A s / Ao < e^ 2 , with 
another pair appearing at A s /Ao < e _7r / 2 . 

There is a direct correspondence between this behav- 
ior of the roots and the dynamical regimes A, B, and 
C observed numerically. The pairing amplitude is sub- 
ject to fast dephasing and tends to zero when L 2 (y) does 
not have complex roots. A pair of complex roots y a = 
±iA a /A s defines the long-time asymptote A(t) w A a . 
Two pairs of roots, y = ±iu\ and y = ±i«2, correspond 
to the parameters A± = (u\ ± U2) A s of the elliptic func- 
tion which defines the asymptotic behavior: 

A(t) =A + dn[A + (t-r ),fe], k = l-A 2 _/A\, (13) 

where the time lag To is a half of the period. As illus- 
trated in Fig0J Ea. H13l) agrees well with A(t) found nu- 
merically. Thus the spectral analysis is in accord with 
the simulation of Bloch dynamics. It confirms the ex- 
istence of the three regimes and also provides the exact 
values A AB = e _7r/2 A and A BC = e 7r / 2 A . 

Finally, we estimate the change of scattering length 
required to cross the A-B and B-C transitions. Us- 
ing the BCS gap in a weakly interacting Fermi gasjj^, 
A = 0.49i?Fe -1 / 9 , g = ^k F \a\, we see that the condi- 
tions A s /A = e ±7r / 2 , written as 1/g — l/g s = ±7r/2, 
translate into 1/a — l/a s = ±kp. At weak coupling 
this corresponds to a small change of scattering length, 
5a/a Ri ±fcj?o, easily achievable for magnetically tunable 
Feshbach resonance. 
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